
#include "auxmath.h"

#include <cassert>

float Math::factorial(const int i) {
   assert(i >= 0);
   if( i == 0 or i == 1) return 1;

   float f = 1;

   for(int j = 2; j < i + 1; ++j) {
      f *= j;
   }

   return f;
}


